One-loop calculations for SUSY processes 



T. Fritzsche and W. HoUik ^ 

'^Max-Planck-Institut fiir Physik (Wcrncr-Heisenberg-Institut) MPP-2004-83 
Fohringer Ring 6, D 80805 Munich, Germany 

Strategy and results for complete one-loop computations in the Minimal Supersymmetric Standard Model are 
reviewed, with applications to the calculation of SUSY mass spectra and SUSY-particle processes. Determination 
of renormalization constants and counterterms are described in the on-shell renormalization scheme, and a trans- 
lation between DR and on-shell parameters is given. As an example, cross sections for chargino and neutralino 
pair production in e^e~ annihilation are presented, complete at the one-loop level. 



1. Introduction 

Experiments at future high-energy colliders 
will be able to discover supersymmetric parti- 
cles and to investigate their properties. A lin- 
ear electron-positron collider will be the best en- 
vironment for precision studies of supersymmet- 
ric models [1], especially of the minimal super- 
symmetric standard model (MSSM). From pre- 
cise measurements of masses and cross sections, 
the fundamental parameters of the MSSM La- 
grangian can be reconstructed [2] to shed light 
on the mechanism of SUSY breaking. Ade- 
quate theoretical predictions matching the exper- 
imental accuracy require the proper inclusion of 
higher-order terms in the calculations of mass 
spectra, cross sections, decay rates 

In this talk we present a brief overview on 
complete one-loop calculations for SUSY pro- 
cesses and supersymmetric mass spectra, where 
all MSSM particles with electroweak couplings 
are included in the virtual states. As a calcula- 
tional frame, the on-shell renormalization scheme 
has been chosen where all particle masses are 
defined as pole masses, i.e. on-shell quantities. 
Cross sections are thus directly related to the 
physical masses of the external particles and of 
the other particles entering the loops. We outline 
the theoretical basis of the calculation and show 
in a few examples the numerical size of the loop 
effects, where we restrict ourselves to the case of 
the CP-conserving MSSM with real parameters. 



The presentation is based on [3,4], for other ap- 
proaches see [5,6]. 

2. Renormalization 

2.1. Charginos and neutralinos 

The bilinear part of the Lagrangian describ- 
ing the chargino/neutralino sector of the MSSM 
involves the ^ parameter, the soft-breaking 
gaugino-mass parameters Mi and M2, and the 
Higgs vacua Vi, which are related to tan/? = 
V2/V1 and to the W mass Mw = gv/2 with 

Renormalization constants are introduced for 
the chargino mass matrix X and for the chargino 



fields X 
X 



1,2) by the transformation 



X + 6X . 



^RX 



(1) 



The matrix 5X consists of the counterterms for 
the parameters in the mass matrix X , 



SX = 



6M2 



V2S{Mw cos/3) 



V2S{Mw sin/3) 



(2) 



The field-renormalization constants SZ^+ and 
are general complex 2x2 matrices. The 
set of renormalization constants in (1) renders 
both S'-matrix elements and Green functions for 
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charginos finite and allows moreover to get the 
matrix of renormalized self energies diagonal on 
the mass shell of each of the chargino mass eigen- 

statcs X]^2- 

The neutralino fields Xi* (* = Ij • ■ • , 4) and the 
mass-matrix Y are - in analogy to the chargino 
case - renormalized by the substitutions 



Y 

Xi 



Y + 5Y , 



(3) 



Besides those parameter counterterms already 
present in (2), the coimterterm matrix 5Y con- 
tains the counterterms for the soft-breaking 
gaugino-mass parameter Mi , for the Z mass and 
the eleetroweak mixing angle, respectively. The 
matrix- valued renormalization constant dZ^jXi is a 
general complex 4x4 matrix. 

Using the on-shcU approach of [3], the pole 
masses of the two charginos, m-+,m^+, and of 
one neutralino, m^o, are considered as input pa- 
rameters, to specify the chargino/neutralino La- 
grangian in terms of physical quantities. This 
is equivalent both to the specification of the pa- 
rameters iJ,,Mi,M2, which are related to the in- 
put masses in the same way as in lowest order, 
as a consequence of the on-shell renormalization 
conditions, and to the definition of the respective 
counterterms. In that way, the tree-level masses 
m-+ and m^o as well as the counterterm matrices 
SX and 5Y are fixed. 

Furthermore, one requires for both, charginos 
as well as neutralinos, that the matrix of the 
renormalized one-particle-irreducible two-point 

^ (2) 

vertex functions T\j becomes diagonal for on- 
shell external momenta. This fixes the non- 
diagonal entries of the field-renormalization ma- 
trices 6Z^+, 6Z\, and SZ^a: their diagonal en- 
tries are determined by normalizing the residues 
of the propagators to be unity. 

2.2. Sfermions 

In the sfcrmion sector, we introduce renormal- 
ization constants for the mass matrices in the L, R 
basis for each sfermion species, , and for the 
mass-eigenstates (s = 1,2 in general, s = L 



for f = v) hy means of the transformation 



r 



{5, 



5Zf 



t) 



(4) 



The matrix is made of the counterterms for 
the parameters in the mass matrix . In the 
case of sneutrino fields, the respective mass ma- 
trix actually consists of only one single element. 
The field-renormalization constants bZ^ are gen- 
eral complex 2x2 matrices, except for the j = v 
case with a single field-renormalization constant 
bZ^ = (5Zl only. 

Applying the renormalization procedure of [4] , 
the sfermion soft-breaking L, i? mass parame- 
ters for the members of each isodoublet, together 
with their renormalization constants, are deter- 
mined via on-shell mass renormalization condi- 
tions for three independent pole masses in the 
case of squarks and two in the case of sleptons. 
The trilinear couplings Af can be related to the 
non-diagonal entries of the field-renormalization 
matrices , and are determined via the re- 
quirement of having diagonal two-point vertex 
functions for on-shell momenta, as in the case 
of charginos and neutralinos. The diagonal en- 
tries of bZ^ are determined by normalizing the 
residues of the propagators to unity. 

2.3. Other sectors 

The formal description of parameter and field 
renormalization in the Standard-Model-like part 
of the MSSM is taken over from [7], yielding for 

the electric charge 

e e = (1 bZ^) e , (5) 

and for the masses of M^, Z, and of the fermions 



bM? 



W,Z I 



rrif mf + bruf. (6) 

The counterterms are determined by the usual on- 
shell conditions in each case; field renormalization 
is treated as described in [7]. 

Renormalization of tan/3 as the ratio of the 
VEVs of the two Higgs fields is done in the DR- 
schcme [81, 



tan/3 : 



V2 



V2_ 
Vi 



5ZhM (7) 
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with the DR-singular parts of the Higgs-field 
renormaUzation constants SZh^ ^- 

3. DR versus OS scheme 

In the DR scheme, divergent 1-loop quantities 

are rcnormalizcd by adding countcrtcrms that are 
proportional to the divergent parts, 

2 

7 + log 47r , 

e 

of the 2- and 3-point vertex functions, regularized 
using the dimensional-reduction method. As a 
consequence, physical observables depend on the 
scale /Udr;. Input variables in the DR scheme are 
a natural choice for GUT-inspired parameter sets 
(e.g. in SUGRA scenarios). 

On the other hand, in the OS scheme the renor- 
malization constants are fixed at physical scales; 
observables arc thus scale independent. The OS 
scheme is convenient for calculations of cross sec- 
tions and decay rates, because masses at Born 
level and in higher order agree (with few excep- 
tions), holding the correct phase-space kinematics 
already in tree- level calculations. 

In order to achieve a translation between DR 
and OS parameters, the following two steps are 
performed, specified here for the quantities fi, M2, 
Ml of the chargino/neutralino sector. 

1. Using n, M2, and Mi in the DR scheme as a 
starting point, the pole masses of three par- 
ticles, e.g. both charginos and the lightest 
neutralino, are calculated at the one- loop 
level. 

2. From those three physical masses the corre- 
sponding parameters in the OS scheme are 
deduced, using tree-level relations (which 
are left unaltered by construction in the OS 
scheme). 

In the example below we show the values (in 
GeV) for the various quantities using the SPSla 
set [9] of input parameters in the DR scheme, 
where shifts up to 10 GeV can occur. 
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DR 


pole 


masses 




OS 




= 352.4 




= 184.3 




= 361.6 


Ml 


= 99.0 




= 387.4 


Ml 


= 102.4 


M2 


= 192.7 




= 99.5 


M2 


= 201.0 



4. e+e production cross sections 

As a concrete example, we discuss the produc- 
tion of chargino and neutralino pairs in e+e^ 
annihilation (see also [10]). For production of 
sfermion pairs we refer to the literature [11]. 

4.1. Born amplitudes 

At lowest order, the amplitude M for chargino- 
pair production can be described by s-channel 
photon and Z-boson exchange and by t-channel 
exchange of a scalar neutrino as displayed in 
the following Feynman diagrams. 

Xi 

/ « \ ^ — 

e 

In the case of neutralino pair production, there is 

no photon exchange at tree level, and i-channel 
exchange is mediated by the two selectrons e*, 
s = 1,2. 




Other diagrams, containing Higgs lines, are al- 
ways negligible. 

4.2. Virtual corrections 

The above set of Born diagrams has to be 
dressed by the corresponding loop contribu- 
tions containing the full particle spectrum of 
the MSSM. The renormalization constants deter- 
mined in section 2 are complete to deliver all 
counterterms required for propagators and ver- 
tices appearing in the amplitudes, and they have 
become part of the FeynArts package for the 
MSSM [12]. The one-loop contributions can be 
classified as follows. 
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• self-energy contributions to the propagators 
for photon, Z boson, and I'e (e'^) sneutrino 
(selectrons); 

• vertex corrections for 7 and Z in the s chan- 
nel and v'f. (e*) in the t channel, 

• box-diagram contributions with double ex- 
change of gauge, Higgs, and SUSY particles 
in the s and t channel. 

4.3. Real photons and "QED corrections" 

Virtual photons attached to external charged 
particles give rise to infrared (IR) divergences in 
the loop diagrams. An IR-finite result is obtained 
by adding rcal-photon brcmsstrahlung integrated 
over the photon phase space. The sum of the 
one-loop contribution to the cross section, cr^"'*, 
and the bremsstrahlung cross section, o-brems^ jg 
IR-finite. Different from standard-fermion pro- 
duction, these photonic contributions cannot be 
removed and treated separately as "QED correc- 
tions", because the one-loop result without the 
virtual photons would not be UV-finite. The 
presence of the photon is required to cancel the di- 
vergence from the photino component of the vir- 
tual neutralinos. 

For cancellation of the IR divergence, it is con- 
venient to split a^'"^"^^ into a (IR-divergent) soft 
part and a (IR-finite) hard part, both depending 
on a soft-photon cutoff AE for the energy of the 
radiated photon. 



brems 



hard / 



(8) 



For practical calculations, in order to avoid 
numerical instabilities, cr'^*''^ is divided into a 
coUinear part, where the photon is within an an- 
gle smaller than A9 with respect to the radiating 
particles, and the complementary non-collinear 
part. 



a 



hard 



= a 



coll 



non— coll 



(9) 



For a^°^^ and (t'^°" analytical results are avail- 
able; cr"°"~'=°'' is calculated numerically with the 
help of DIVOMME and CUHRE, multidimensional nu- 
merical integration routines that are part of the 
CUBA-hbrary [13]. 

As already pointed out above, there is no dia- 
grammatic way to disentangle QED-like photonic 



virtual contributions from the MSSM-specific 
parts. One can, however, isolate the universal 
and leading QED terms in a™^ + (7^°^* resulting 
from photons collinear to the incoming e^, which 

contain the large logarithm = log • 
The separation 



virt I -.soft I ^ 

\ (J — C/ \ (^remainder : 



(10) 



a 

a = — 

TT 



(ie-l)log 



4AE 



s 2 ' 



^Born 



identifies a one-loop contribution cTremainder that 
is IR finite and free of large universal QED terms, 
"■remainder Contains the MSSM-specific radiative 
corrections, whereas the subtracted part a in (10) 
together with the hard bremsstrahlung part from 
initial-state radiation can be considered as a con- 
tribution of the type "QED corrections" , 



"■QED 



hard 
'^ISR 



(11) 



It is independent of the auxiliary cut AE, and 
it contains all large logarithms from virtual, soft, 
and hard photons. 

With this reordering, the complete cross sec- 
tion at the one-loop level can be written as fol- 
lows, 

^1— loop „Born ■ ^virt i ^soft i ^hard /i o^ 

a ^ = a +a +a +a (izj 

= (7^°''" -|- (TQED + ^remainder + Acr. 

The small piece Aa accounts for final-state ra- 
diation in case of chargino production, and for 
photon radiation from the f-channel selectrons in 
case of neutralino production. 

4.4. Results 

For illustration, we choose the example of the 
SPSla parameter set [9] listed partially in sec- 
tion 3. The DR parameters arc first converted 
into on-shell quantities and pole masses, from 
which the cross sections are calculated, both at 
trcx' \v\iA and at higher order. The following 
figures contain the integrated cross sections for 
XiXi and for X1X2 production in e+e" annihi- 
lation for unpolarized beams, as functions of the 
CMS energy. Besides Born cross sections and full 
one-loop results, also the contributions (Tremainder 
and (TQED are shown separately. 
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